function zhjsuan2301
x=linspace(0,2*pi,50);
y=f1(x);
plot(x,y,'-or','LineWidth',2)
hold on
plot(x,f2(x),'s','LineWidth',1,'Markerface','b','Markeredge','r')
ezplot('abs(x).^(4/5)+abs(y).^(4/5)-1',[-1 1])
hold off
print -djpeg90 ͼƬ1

function y=f2(x)
y=f1(x).^2;
function y=f1(x)
y=sin(x).*cos(x);